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Quantum information processing has been efFectively demonstrated on a small number of qubits 
by nuclear magnètic resonance. An important subroutine in any computing is the readout of the out- 
put. "Spectral implementation" originally suggested by Z.L. Madi, R. Bruschweiler and R.R. Ernst, 
[J. Chem. Phys. 109, 10603 (1999)], provides an elegant method of readout with the use of an extra 
'observer' qubit. At the end of computation, detection of the observer qubit provides the output via 
the multiplet structure of its spectrum. In "spectral implementation" by two-dimensional experi- 
ment the observer qubit retains the memory of input state during computation, thereby providing 
correlated information on input and output, in the same spectrum. "Spectral implementation" of 
Grover's search algorithm, approximate quantum counting, a modified version of Berstein-Vazirani 
problem, and Hogg's algorithm is demonstrated here in three and four-qubit systems. 

I. INTRODUCTION 



' In 1982 Feynmann pointed out that it would be more efíicient to simulate the behavior of a 

^ .quantum system using a quantum, rather than a classical device [1]. The idea of a purpose-built 
FT' 

^ quantum computer which could simulate the physical behavior of a quantum system as well as 
• i— i , 

^ 'perform certain tasks much faster than classical computer, attracted immediate attention [2,3]. The 

h ; 

5t .theory of such quantum comp úters is now well understood and several quantum algorithms like 
Deutsch-Jozsa (DJ) algorithm, Grover's search algorithm, Shor's factorization algorithm, Berstein- 
Vazirani problem, Hogg's algorithm and quantum counting have been developed [4-10]. 

However, building a realistic large scale quantum computer has been extremely challenging [11,12]. 
Various devices are being examined for building a quantum information processing (QIP) device 
which is coherent and unitary [11]. Among these, nuclear magnètic resonance (NMR) has shown 
great promise by demonstrating several quantum algorithms and other QIP tasks on small-scale 
devices. [13-28]. The last step in any quantum information processing task is the "readout" of 
the output. Typically in NMR, the readout is obtained by selectively detecting spins [29], or by 
mapping out the full density matrix [30-32]. It was first pointed out by Ernst et. al. [33] that it is 
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advantageous from the spectroscopic viewpoint that quantum states can be assigned to individual 
spectral lines, corresponding to transitions between energy levels rather than to the energy levels 
themselves [33]. However, for such an advantage one has to use an extra qubit called "observer" qubit. 
After computation the readout is obtained by detecting only the observer qubit, whose multiplet 
spectrum provides the result of the computation carried out on the work qubits. Such a "spectral 
implementation" of a quantum computer was demonstrated by implementation of some lògic gates by 
one- and two-dimensional NMR [33]. Later, "spectral implementation" of a complete set of lògic gates 
and DJ-algorithm [34], Berstein-Vazirani problem [35] and quantum Fourier transform [36] has also 
been implemented by NMR. In this work we extend this range by spectrally implementing Grover's 
search algorithm, approximate quantum counting, a modified version of Berstein-Vazirani problem, 
and Hogg's algorithm. All the algorithms are implemented by both one- and two-dimensional NMR. 
To the best of our knowledge this is the first "spectral implementation" of these algorithms. 



A convenient representation of the density matrices of pure states in Liouville space can be obtained 
by the polarization operators for each qubit (k) [33,37], 
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For example, the density matrix of a pure state 1 00) + |11) can be expressed as 
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|00)(00| + |11)(11| + |00)(11| + |11)(00| 
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The scheme of "spectral implementation" of one-dimensional (1D) and two-dimensional (2D) NMR 
is respectively given in figure l(a) and l(b). We start with the thermal equilibrium density matrix 
7° + 1\ + 1% + .. + 1^ and in the preparation period we create density matrix of the form , 
where = (I^I^.J^ - I^I^.J^ )/2. In this state the last N-l qubits are simultaneously 

in pseudopure state (PPS) [13] in two distinct domains of energy levels, in which the observer qubit 
is in state |0) and |1) respectively. Such a state is known as sub-system pseudopure state [14]. This 
is further elaborated in figure 2. 

The schematic diagram of the energy levels and the spectrum of the observer qubit in a three 
qubit system, where the first qubit is the observer qubit, is given in figure 2. Figure 2(a) shows 
the equilibrium deviation populations (populations in excess of uniform background population) of 
various energy levels and figure 2(b), the equilibrium spectrum of the observer qubit obtained after a 
(n/2) detection pulse. Each of the spectral lines in the multiplet, correspond to the state of the other 
qubits. The energy level diagram along with the deviation populations after creating the desired 
initial state of I^Iq^o * s given in 2(c). The corresponding observer qubit spectrum has a single line, 
that of 1 00) , indicating that the other qubits are in 1 00) state. 

Typically after computation, the density matrix is of the form 7°Jq; 1 Jq/ 1 ..J(w 1 , where the subscript 
0/1 means that the particular qubit is either in or 1 state. A subsequent (vr/2)° pulse on the 
observer (7°) qubit creates single quantum coherences of the form 7°7o/i-^o/i··-^o/i' wri ich gives a 
single line in the spectrum corresponding to the output state of other qubits. An example for the 
3-qubit system is given in figure 2. Let us assume that we start with the initial 1 00) pseudopure 
state of the qubits (other than observer qubit) and after some computation let the output state be 
|11). After such a computation, the deviation populations and spectrum of observer qubit are given 
respectively in figure 2(e) and 2(f). 

In some algorithms however, the output is a superposition of múltiple states. Then, the output 
density matrix will have non-zero populations in all the output states and the coherences between 
them. The spectrum of the observer qubits will thus have múltiple lines, corresponding to all the 
output states. For example, in the 3-qubit system, if the output state of the work qubits is 1 00) + 1 11), 
the density matrix is of the form J°#ig + I° z l\ll + 7°/+/+ + 7 2 °/Í/ 2 . After the (tt/2)° detection 
pulse on the observer qubit the single quantum coherences of the terms I^IqIq and I°J\ll will be 
detected. The spectrum of the observer qubit will show two lines corresponding to the states of 1 00) 
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and |11) of the other qubits. The coherences will be converted into múltiple quantum coherences 
which are not detected directly in NMR. Hence, the "spectral implementation" gives a measure of 
the deviation populations or probabilities of each state but does not measure the coherences, which 
if required can be measured by state tomography [30-32]. 

A two-dimensional experiment for "spectral implementation" provides the input and output in the 
same spectrum. The pulse sequence for the two-dimensional experiment of "spectral implementa- 
tion" is given in figure l(b). Suppose a computation starts with the input of 1 00. .0) and end with an 
output of 1 11.. 1) state. After preparation of the initial I®IqIq..Iq state the application of the pulse 
sequence of figure l(b) can be analyzed in the following steps: 

r0rlr2 tN (*72)g T° T 1 T 2 T N 
I z I I 0·· I I x I I 0-· I 

W4«- Í2 ^íi?··/f r ««(ü;a ) ..o*i)oo a (a;? 1 .. 1 í 2 ), (4) 

where Uq and co^ A are respectively the freqüències of the |00..0) and |ll..l) transitions of the 
observer qubit 1°, (7r/2)° is a (ir/2) rotation of the observer qubit (7°) about y-axis, G z is the gradient 
pulse and Comp is the computation performed on the work qubits. It may be noted that the signal 
from the observer qubit is modulated by the freqüències corresponding to both the input and the 
output states of the work qubits. A series of experiments are performed with systematic increment of 
the íi period followed by detection of the observer qubit's signal. The collected two-dimensional time 
domain data set s(íi,í 2 ) is double Fourier transformed yielding a two-dimensional frequency domain 
spectrum S(lüi, cj 2 ), which contains along o>i the input states of work qubits before computation and 
along üü 2 , the output state of work qubits after computation. 



III. GROVER'S SEARCH ALGORITHM 



Grover's search algorithm can search an unsorted database of size N in 0(y/~N) steps while a classi- 
cal search would require O(N) steps [5]. Grover's search algorithm has been earlier demonstrated by 
NMR [18,19]. The quantum circuit for implementing Grover's search algorithm on two qubit system 
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is given in figure 3(a). The algorithm starts from a 1 00) pseudopure state. A uniform superposition 
of all states are created by the initial Hadamard gates (H). Then the sign of the searched state "x" 
is inverted by the oracle through the operator 

U x = I-2\x)(x\. (5) 

An inversion about mean is performed on all the states by a diffusion operator HUqqH , where 

f/ o = /-2|00)(00|. (6) 

For an N-sized database the algorithm requires 0(y/N) iterations of U x HUqqH. For a 2-qubit system 
with four states, only one iteration is required. We have implemented this algorithm on the two 
qubits of a three qubit system with the third qubit acting as the observer qubit. The three qubit 
system chosen for this purpose is 4-fluro 7-nitro benzofuran (dissolved in CDCI3), which comprises 
of a two protons ( 1 H) and a flourine ( 19 F). The chemical structure of the molecule along with the 
equilibrium proton and fluorine spectrum is given in figure 4(a). We have chosen the fluorine spin 
as the observer qubit. The Hamiltonian of the system is 

n = jz2^ t r z + Y.^Wl (7) 

i=0 i>j 

where z/j are the resonance freqüències of various spins and are the indirect couplings. The 
experiments were performed at a field of 11.4 Tesla in a Bruker DRX500 spectrometer. At the 
magnètic field of 11.4 Tesla, the resonant frequency of proton is 500.13 MHz and that of fluorine is 
470.59 MHz. The frequency difference between the two protons is 646 Hz. The J-couplings are J01— 
-3.84 Hz, J 2=8.01 Hz and J\i— 8.07 Hz. The l· l·l transmitter frequency is set at the center of the 
proton spectrum. 

The required initial state of J°|00)(00| was prepared by the method of pair of pseudopure states 
(POPS), originally suggested by Fung [38,39]. The method requires two population distributions, (i) 
equilibrium populations and (ii) population distribution after a selective (7r)-pulse on 1 000) <-> 1 100) . 
Subtraction of (ii) from (i) effectively gives the initial state of 7°|00)(00| (figure 4(b) corresponding 
to the schematic PPS of figure 2(c)). It might be noted that the method of creation of sub-system 
pseudopure states from cat-states can also be adopted for creation of this initial state [40]. 

The Hadamard gates are implemented by (7r/2)^(7r) : J,' 2 pulse (pulses are applied from left 
to right) [18], where (9) x ' 2 denotes a 0-angle pulse (rotation) on l st and 2 níi -qubit about the 
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rr-axis. The Uqq operator is a controlled phase gate which can be implemented by the se- 
quence [(r/2)(7r)i' 2 (r/2)(7r)i' 2 ] [(tt/2) ^ (tt/2) Í^(tt/2) J' 2 ] , where r = 1/2 J 12 [18]. The sequence 
[(t/2)(tt)]: 2 (t/2)('k)1: 2 ] evolves the system only under the Ji2-coupling and refocuses all other cou- 
plings and proton chemical shifts [37], whereas the [(7i/2) 1 _ly(7r/2) 1 _l 2 : (7r/2)y' 2 ] is a composite z-rotation 
on both the qubits [41]. Similarly, the other phase gates can be constructed as [18], 

U 01 = [(r/2)(7r)^(r/2)( 7 r)^][(7r/2)^(7r/2)i(7r/2) 2 _ :c ( 7 r/2);' 2 ] 
U W = [(r/2)( 7 r)^ 2 (r/2)(7r)^ 2 ][( 7 r/2)^(V2)^( 7 r/2) 2 (7r/2);' 2 ] 
U n = [(r/2)( 7 r)^ 2 (r/2)( 7 r)^ 2 ][( 7 r/2)^( 7 r/2)^ 2 ( 7 r/2);' 2 ] (8) 

The pulses which are simultaneously applied on both the qubits are achieved by hard pulses. How- 
ever, some gates require selective excitation of qubits. Since the resonance freqüències of the two 
protons are relatively close to each other, selective excitation of a particular proton qubit requires 
long low-power pulses, which introduce significant errors in the computation [17,19]. Fortunately, 
in case there are two homonuclear qubits, the selective pulses can be substituted by hard pulses 
and delays using the variation of "jump-and-return" sequence [42], as demonstrated by Jones et. al. 
[43]. For example, the pulse sequence of Uoi gate requires (t: /2) 1 x (t: /2) 2 _ x at one point. This can be 
achieved by using the identity [41] 

(n/2)_ y (n/2) ±z (n/2) y = (n/2) ±x . (9) 

If the proton transmitter frequency is set at the center of the spectrum, then v\ — — u 2 — v, and, a 
delay of (l/4i/) evolves the two protons under the Zeeman Hamiltonian of 2txv(I\ — I 2 ) to give the 
intermediate (vr/2) ±2 rotation of Eq. [9] . Hence, 

(^2)^/2)^ = (n/2U - (1/41/) - (tt/2),. (10) 

Similarly, the pulse (7r/2) 1 i ;£ .(7r/2) 2 required for U\o gate, can be achieved by 

(7r/2)^(vr/2) 2 = (tt/2), - (l/4i/) - (n/2)_ y (11) 

In principle however, the evolution under J-coupling during (l/4i/) would lead to some non-ideal 
characteristics [44], which is minimal in our system, since the ratio of maximum J-coupling to 
chemical shift frequency difference ~ 1:80. This error is significantly less than the error introduced 
due to evolution under interna! Hamiltonian during low power long duration qubit selective pulses. 
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After application of the quantum circuit in figure 3(a) on the initial state of Iq 1 00) (00 1 , the observer 
qubit was detected by a (tt/2) pulse. From the obtained spectrum given in figure 5(a), one can identify 
the searched state (\x)) directly. The two-dimensional experiment for "spectral implementation" 
has the added advantage that the input and output can be identified in a single spectrum. The 2D 
experiment of figure l(b) was carried out, where during the computation period, the quantum circuit 
of figure 3, was implemented on li and I 2 . The resultant spectrum given in figure 5(b), shows the 
input and output in each case. For example, when \x) = |11), a cross-peak at the frequency of 1 00) 
transition along lú\ to that of the |11) transition along a; 2 , identifies the input as 1 00) and the output 
as |11). The 2D spectra in figure 5(b) contains the initial state of 1 00) and the searched state of 
1 00) , |01), 1 10) and |11). In the 1D spectrum of Fig. 5(a), the intial PPS state has to be ascertained 
independently prior to the implementation of search algorithm. 



The search problem may be thought as finding fc-entries out of N, which satisfy the condition 
f(x) = 1. For the other, N — k entries, f(x) = 0. While Grover's search algorithm searches 
these k-items (one at a time), quantum counting finds out the value of k [9,10]. This has extreme 
importance because in case of múltiple solutions, the required number of Grover's iteration scales as 



0(yN/k) [12]. Hence finding out the number of existing solutions speeds up the search procedure. 
Moreover, the fact that counting can find out whether the number of solutions is zero or finite, 
makes it applicable to the non-deterministic (NP)-complete search problems, where it is important 
to know whether solution exists for a given search problem [12]. Approximate quantum counting 
has been demonstrated using NMR by Jones and Mosca [43]. In this work we provide a "spectral 
implementation" of approximate quantum counting in the three qubit system of 4-fluro 7-nitro 
benzofuran, where the I 2 is the target qubit, I 1 the control qubit and the 1°, the observed qubit. 

The working of counting algorithm, as detailed by Jones and Mosca [43], is as follows. Count- 
ing algorithm can be thought of as a method for estimating the eigenvalue of Grover's iteration 
G = HUoH^Uf, where U = I - 2|00..0) (00..0| and U f transforms \x) to (-iyW +1 \x). Starting 
from the initial 1 00. .0) state, an initial Hadamard on target qubit creates an uniform superposition 
H\00..0) = (\íp + ) + \i>-))/\/2, where and are two eigenvectors of G [43]. These two eigen- 
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vectors are with eigenvalues of e ±t<t>k , where sin((f>k/2) = yk/N. An uniform superposition of the 
control qubit is also created. The application of controlled G produces the result 

K) = (|0)+e^|l))|^ + )/ V / (2). (12) 

If r iterations are performed, then the state is 

|^) = (|0)+e^|l))l^ + )/ v / (2). (13) 

A second Hadamard gate on the control qubit produces 

\V Í+ ) = [(1 + e**")\0) + (1 - e^)|l)]|<M/2. (14) 

A similar result will happen in the case of \tp-). At the end, the final state \ipf) will be an entangled 
state of the control and target qubits, except when k = or k = N [9,10,43]. 

Jones et. al. have implemented the quantum circuit of figure 3(b) in a two-qubit system, measured 
the signal from control qubit, thereby tracing the target qubit, and shown that the signal assumes 
a sinusoidal behavior with r whose frequency depend on 0& [43]. We have instead, started from the 
initial J°|00)(00| state and inferred the result of counting from the spectrum of observer qubit. For 
a two qubit case only one Grover's iteration is sumcient to get the result [43] . Given in Table 1 are 
the count k, their corresponding 0, the Uf operators and final state of the system for r — 1. Note 
that for k — 0, the final state is \ipf) = 1 00) and for k = N = 2 the final state is \ipf) = |10). For 
k — 1 the output states are in entangled form of all the states, 1 00) , 1 01) , 1 10) and |11). 

Starting with the initial state of J° 1 00) (00 1, we implemented the quantum circuit of figure 4(b). The 
controlled controlled U and controlled U f 01 have the same operator as that of two-qubit controlled 
phase gate U w implemented in Grover's algorithm (section II), whose corresponding pulse sequence 
is given in Eq. [11]. Uf 10 has the same operator and pulse sequence as that of Un in Eq. [11]. C// n 
is an identity operator and required no pulses. Uf 00 required a rotation. This (tt)1 rotation was 
implemented with hard pulses and evolution under Zeeman Hamiltonian, 

= (V4^) - [(7r/2)^(7r/2)i' 2 (7r/2);. 2 ]. (15) 

During the delay (l/4i/) the system evolves under the Zeeman Hamiltonian to acquire a rota- 
tion of (tc/2)1(tc/2) 2 _ z . The subsequent composite z-pulse was applied on both qubits, (tx/2) 1 ^ 2 = 
[(7r/2)^(7r/2)^ 2 (7r/2)^' 2 ], which cancels the rotation of second qubit but adds to the rotation of first 
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qubit to give an effective (n/2)\ rotation. It may be noticed that there are two pseudo-Hadamard 
gates on second qubit which require spin-selective pulses since h=(n/2) y and hr 1 =(n /2)_ y . However, 
these pulses can also be performed by hard pulses and evolution under Zeeman Hamiltonian using 
the "jump-and-return" lògic [43]. 

(tt/2); = (n/2) 1 / - (I/81/) - (7r/2)^(7r/4)J' 2 , 

(n/2)*_ y = (n/2) 1 / - (l/8u) - (n /2)/(n / A) 1 /. (16) 

After implementing the quantum circuit of figure 3(b), the observer qubit was measured. The 
observer qubit's spectrum given in figure 6(a), shows four lines for k = 1 (/oi and fio). For k — (foo), 
only 1 00) transition and for k = 2 (fn), only 1 10) transition is observed. The 2D-spectrum of figure 
6(b) contains correlation of the output state with the initial 1 00) pseudopure state, and confirms the 
same result. 

V. BERSTEIN-VAZIRANI PROBLEM 

Berstein and Vazirani considered the problem of determining a n-bit string "a" [8]. Classically 
each query would yield one bit of information and hence would require n-queries to the database. 
However, Berstein and Vazirani showed that a quantum algorithm can solve the problem with one 
quantum query [8]. For this purpose, the oracle has to compute a function f a (x) = a.x . The scheme 
proposed by Berstein and Vazirani required an ancillary qubit and determined a n-qubit string with 
n+1 qubits, which has been demonstrated by NMR recently [36]. However, Du and his co-workers 
had simplified the scheme such that the ancillary qubit was not required [46]. We have implemented 
the Du-scheme, since it has the advantage of determining a n-qubit string with n-qubit system. 
The quantum circuit of a two-qubit implementation is given in figure 3(c). Starting from |0) n , the 
Hadamard gates create an uniform superposition 

The U a operator transforms \x) — > (—l)^ a ^\x). The unitary operator U a can be decomposed into 
direct products of single-qubit operations [46] 

u a = u 1 ®u 2 <g> ...®u n , 
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XP = I, ai = 
= cr z , = 1 

Operation of C/ a creates a new state |t/> 2 ) of the form, 

w = w = o^E(-im (19) 

The final state after the subsequent Hadamard operation is 

2 n_ 12 n_ 1 

l^> = W 2 ) = ^ E E (-i) M (-inv> (20) 

However, since Ef^-l)^-!) 3 ^ = ^ [2], |^ 3 ) = |a) [8,43]. 

The algorithm was implemented to determine a two-qubit string by "spectral implementation" 
using three qubit system of 4-fluro 7-nitro benzofuran. After creating POPS, followed by Hadamard 
pulses, the operator U a was applied for \a) = 1 00) , |01), 1 10} and |11). C/oo is unity operator and 
does not require any pulse. Uio is al, which requires a (ir)l rotation. Once again, the (ir)l rotation 
was implemented using the pulse sequence of Eq. [15]. Similarly, C/ i was implemented by 

= WÍ' 2 (1/^)WÍ' 2 - [(7r/2)^(7r/2)^(7r/2)f]. (21) 

Un is which can be achieved by a composite z-pulse of (7r)?;' 2 = (7r/2) 1 L^(7r)J.' 2 (7r/2)^ 2 . After 

application of the final Hadamard pulses, the observer qubit was detection by a (ir/2) pulse. The 
obtained spectrum given in figure 7(a), clearly determines the 2-bit string in each case. The result 
of 2D experiment is given in figure 7(b). The 2D spectrum correlates the input 1 00) to the output 
in each case. 

The above algorithm was also implemented to determine a three-qubit string by "spectral im- 
plementation" using a four-qubit system. The molecule 2-3 diflouro 6-nitrophenol (dissolved in 
CDCI3+I drop D 2 0) has 4 weakly coupled spin-1/2 nuclei. The proton of the phenol group is 
exchanged with the D 2 0. The two remaining protons and the two fluorine nuclei constitute the 
four-qubit system. The equilibrium spectrum of each nucleus is given in Fig. 8 (a). In a 500 MHz 
NMR spectrometer, the chemical shift difference between the two Fluorine spins is 16 kHz while 
that between the two protons is 560 Hz. The couplings range from 19.13 Hz to -2.4 Hz. 
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The operators and pulse sequences required for each string of a three qubit system is given in Table 
2. Since the chemical shift difference between the two fluorine spins are considerably large (16kHz), 
selective pulses do not introduce significant errors. The pulses on fluorine spin J 3 were achieved by 
Gaussian shaped selective pulses of 12.5/xs duration. The proton transmitter frequency was kept at 
the middle of the spectrum and the selective z-pulses on protons were applied in similar lògic as 
in the two-qubit case (Eq. [15] and [21]). Hard pulses were applied when both protons had to be 
pulsed simultaneously. The algorithm was implemented starting from the initial state of 7°|000) (000 1 
(figure 8(b)) and finally the observed qubit was measured by selective Gaussian shaped (ir/2)y pulse. 
The only transition present in each spectrum given in figure 9(a) indicates the corresponding string. 
The 2D experimental spectra given in figure 9(b) verify the same results, correlating the input state 
of [000) in each case. 

VI. HOGG'S ALGORITHM 

Satisfiability (SAT) problem is one of the nondeterministic polynomial (NP) combinatorial search 
problems [7]. SAT problem consists of a logical formula in n variables, Vi,V 2 , ...,V n [7]. One has to 
find an assigment (true or false) for each variable Vi, such that it makes the formula true. The logical 
formula can be expressed in various equivalent forms, as conjunction of clauses, where a clause in 
a disjunction of some variables. A clause with k variables is false for exactly one set of vàlues for 
its variables but true for the other 2 k — 1 sets. An example of clause for k=3 is V\ OR Vi OR V3, 
where the clause is false for only V\ — V2 — V3 = false. Only the assignments which satisfy all the 
clauses are considered as solutions [7]. 

While the number of steps required by a classical algorithm increase linearly with the size of 
the variables [7], Hogg's algorithm can solve 1-SAT and maximally constrained fc-SAT problems 
in a single step, whatever be the size of the variable [7]. Hogg's algorithm has been succesfully 
implemented by NMR in a three qubit system [45]. Here we demonstrate Hogg's algorithm by 
"spectral implementation" . The Hogg's algorithm starts by creating a uniform superposition of 
states by initializing from |^ ) = \0)® n and applying Hadamard gate on all qubits = H® n \^) ) = 
2~ n í 2 Yj S \ s ) ■ Let m be the number of clauses in the 1-SAT logical formula. Then the unitary 
operations UR are consecutively applied to yield the final state \tpf) = UR\ipi) where U and R 
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defined as follows. R adjusts the phases of \s) depending on the conflicts c of the different assignments 
in the superposition of s, ranging from to m. R is a diagonal matrix of the form 

R ss = V2cos[(2c - 1)tt/4], for even m 

= i c for odd m (22) 

The operator U mixes the amplitudes from different assignments with elements U rs , depending the 
Hamming distance d between r and s. U is of the form 

— 2~ ( "~ 1)/2 cos[(n — m + 1 — 2d)n/4\ for even m 

= 2 -n/2^(n-m)/^_ i y for m ^3) 

where d — \r\ + \s\ — \r A s\ is the Hamming distance between r and s, i.e. number of positions at 
which their vàlues differ. U can be decomposed into HYH where H is the Hadamard matrix and Y 
is a diagonal matrix of the form 

r r r = 7(7) = lh = V2cos[(m — 2h — l)7r/4] for even m 

= i h e -inm/4 for odd m ^ 

where h — \r\, and hence r rr depend on the number of 1-bits in each assignment. For a detailed 
description of the working of the algorithm see ref [7]. Hence the Hogg's quantum starts with the 
initial 1 00) state and reaches the desired output state \ipf) by 

= URH \00) (25) 

and the final step is to measure the output state \tpf). We have observed the final state through 
detection of the observer qubit. 

The quantum circuit of Hogg's algorithm for a three-qubit system is given in figure 3(d). While 
implementing the corresponding pulse sequence, the consecutive pulses of opposite phases cancel out, 
yielding a simplified sequence [45]. The m — 1 and 3 clauses, their lògic formulae and the reduced 
pulse sequences are given in table 3 [45]. Only the m—1 and 3 cases are demonstrated here for ease 
of implementation. The selective pulses on fluorine spin were achieved using Gaussian shaped pulses. 
In protons, the selective pulses were achieved by hard pulses concatenated with Zeeman evolution 
(Eq. [10], [11] and [16]). For example, while implementing Vi, 
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Vi : [(7r)i][(7r/2);][(7r/2);] = [(tt/2)^ - l/4i/ - (^2)^(^/2)^] 

[(tt/2)^ - I/81/ - (7r/2)^( 7 r/4)f][( 7 r/2)J] 
= [(tt/2)^ - 1/41/ - (tt/2)^ 2 - I/81/ - (7r/2)^(7r/4)J' 2 ][(7r/2)J] (26) 

The spectra obtained by the one-dimensional experiment is given in figure 10 while the spectra 
obtained in 2D-experiment is given in figure 11. The spectra of observer qubit clearly identifies the 
desired outputs of table 3. For example, in the case of V 1: the output has all the states that satisfy 
the condition that l st qubit is |1) or false; namely |001), 1 01 1 ) , 1 101) and 1 1 1 1 ) (read the order of 
qubits from right to left). Similarly, for V 3 A V 2 A V 1 , the output consists of the sole answer 1 111) , 
which satisfies the condition that all the qubits are in state |1). 

VII. CONCLUSION 

We have demonstrated "spectral implementation" of several quantum algorithms by one- and 
two-dimensional NMR. Provided future quantum computers run with high fidelity, "spectral imple- 
mentation" delivers an aphoristic method of readout. Though it requires the use of an observer 
qubit, this qubit also helps in creating a pseudopure state by non-scalable and effective methods 
like POPS [38]. With the essentiality that the observer qubit has resolved spectrum, the principle 
of "spectral implementation" is applicable to higher qubit systems without increasing complexity. 
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Table 1 :The various possible the count k for a two-qubit system, their corresponding 0, the Uf 
operators and final state of the system \ip ou t P ut)- 
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Table 2: The operators and pulse sequences required for determination of each string in a three 
qubit system 



string 


Operator 


Pulse sequence 


|000) 
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no pulse 
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o* 
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|011) 
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1 100) 
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1 101) 
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M'M 
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Table 3: Logic formulae for m=l or m=3 in a 3-qubit system, corresponding pulse sequences, and 
theoritical solutions [45]. Read the order of qubits from right to left. 
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Figure Captions 

Figure 1. Experimental protocol for "spectral implementation" of quantum algorothms [33]. (a) 
One-dimensional experiment. The first stage is to create an subsystem pseudopure state of the 
type I®IqIq..Iq , followed by computation on I l ..I N qubits. Finally the transitions of the observer 
qubit 1° are detected by a 90° pulse. (b) Two-dimensional experiment. After the creation of initial 
I^IqIq.Jq subsystem PPS, the observer qubit is flipped by 90° pulse to transverse magnetization 
and allowed to evolve for a time t\. During ti, the transitions of the observer qubit modulate with 
the freqüències characterized by the input state of the other N qubits. A subsequent 90°_ y brings the 
magnetization back to longitudinal direction. The computation is performed on the 7 1 ..7 JV qubits. 
The transitions of the observer qubit are finally detected by a 90° pulse. A series of experiments 
are performed with systematic increment of the t± period and the collected 2D data set s(ti,t 2 ) is 
Fourier transformed to get the 2D spectrum S(ui,U2). 

Figure 2. Schematic diagram of the energy levels and the spectrum of the observer qubit at 
different stages of "spectral implementation". (a) Deviation equlibrium populations. The dotted 
arrows denote the transitions of observer qubit. (b) Equilibrium spectrum of observer qubit shown 
by stick diagram. Each transition of the spectrum correspond to the state of other qubits, which 
are given above each line. (c) Deviation populations after creation of 1 00) subsystem pseudopure 
state by POPS. Populations of only |00) eigenstate is non-zero in the two distinct domains of energy 
levels, where observer qubit is respectively in state |0) and |1). (d) Spectrum of observer qubit 
after creation of POPS. The dots denote null intensity. (e) Deviation populations after a typical 
computation whose output is |11). (f) Spectrum of observer qubit after such a computation. 

Figure 3. The quantum circuits of various algorithms. (a) Quantum circuit for implementa- 
tion of Grover's search algorithm in a 2-qubit system. (b) Quantum circuit for implementation 
of approximate quantum counting in a 2-qubit system. (c) Quantum circuit for implementation of 
Berstein-Vazirani problem. (d) Quantum circuit for implementation of Hogg's algorithm in a 3-qubit 
system. 

Figure 4. (a) Chemical structure and equilibrium spectrum of 4-fluro 7-nitro benzofuran. The 
J-coupling vàlues are Joi= -3.84 Hz, J 2=8.01 Hz and Jyi= 8.07 Hz. The peak denoted by asterisk 
(*) belongs to solvent, (b) Spectra after creation of POPS. A Gaussian shaped selective pulse of 
500ms duration was applied on the |000) <-> 1 100) transition and the resultant spectra is subtracted 
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from the equlibrium spectra of figure (a) to yield (b). [See figure 2(c)] 

Figure 5. (a) "Spectral implementation" of Grover's search algorithm by 1D experiment. After 
computation, the observer qubit is detected by a non-selective pulse of 14 fis. 4x1024 data points 
were collected and zero-filled to 8x1024 before Fourier transform. The observer qubit's spectra 
shows only the transition corresponding to the searched state (\x)) with non-zero intensity. (b) 
"spectral implementation" of Grover's search algorithm by 2D experiment. A 2D data set of 256 x 16 
(í 2 x ti)was collected and zero-filled to 1024x256. It may be noticed that the total size of the raw 
2D dataset is of the same size as that of the 1D experiment. The doubly Fourier transformed spectra 
gives the input state along uj\ and output state along uj 2 . 

Figure 6. (a) "Spectral implementation" of approximate quantum counting by 1D experiment. 
4x1024 data points were collected and zero-filled to 8x1024 before Fourier transform. The observer 
qubit's spectra show the transitions corresponding to the ouput state. Hence the various cases of 
k=0 (/oo), k=l (/oi and / i), and k=2 (/n) can be easily identified from the spectra. (b) "spectral 
implementation" of approximate quantum counting by 2D experiment. A 2D data set of 256x16 
(t 2 x ti) was collected and zero-filled to 1024x256. The Fourier transformed spectra gives the output 
state as well as the input state. 

Figure 7. (a) "Spectral implementation" of Berstein-Vazirani problem in a 2-qubit system. The 
observer qubit's spectra shows the transitions corresponding to the bit string. The strings a = 
00, a — 01, a — 10 and a — 11 can be identified directly from the spectra. (b) 2D "spectral 
implementation" of Berstein-Vazirani problem. A 2D data set of 256x16 (t 2 x í^was collected and 
zero-filled to 1024x256 before Fourier transform. The Fourier transformed spectra gives the bit 
string against the input state in each case. 

Figure 8. (a) Chemical structure and equilibrium spectrum of 4-fluro 7-nitro benzofuran. The 
J-coupling vàlues are J 01 — 5.23 Hz, J 02 = 8.85 Hz, J 03 = 19.1 Hz, J 12 — 9.76 Hz, J 13 — -2.4 Hz and 
J 2 3= 6.81 Hz. (b) Spectra after creation of POPS. A Gaussian shaped selective pulse of 500ms 
duration was applied on the |0000) <-> 1 1000) transition and the resultant spectra is subtracted from 
the equlibrium spectra of figure (a) to yield (b). 

Figure 9. (a) "Spectral implementation" of Berstein-Vazirani problem in a 3-qubit system. After 
computation, the observer qubit is detected by a selective pulse of 12.5us duration. The observer 
qubit's spectra show the transitions corresponding to the bit string. The eight possible strings of a — 
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000, a = 001. ..a =111 can be identified directly from the spectra. (b) 2D "spectral implementation" 
of Berstein-Vazirani problem in the three-qubit case. A 2D data set of 256x24 (í 2 x t{) was collected 
and zero-filled to 1024x256 before Fourier transform. The Fourier transformed spectra give the 
various bit strings along with the input state in each case. 

Figure 10. One-dimensional "spectral implementation" of Hogg's algorithm in a 3-qubit system. 
After computation, the observer qubit is detected by a selective pulse of 12.5 fis duration. The 
observer qubits spectra clearly show the output states corresponding to various logical formulae of 
table 3. (a) contains the spectra corresponding to m=l and (b) to m=3. 

Figure 11. Two-dimensional "spectral implementation" of Hogg's algorithm in a 3-qubit system. 
The two-dimensional spectra provides the output states corresponding to various logical formulae of 
table 3. (a) contains the spectra corresponding to m=l and (b) to m=3. 
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(d) Hogg' s algorithm 
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